Theory of the Three Dimensional Quantum Hall Effect in Graphite 
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We predict the existence of a three dimensional quantum Hall effect plateau in a graphite crystal 
subject to a magnetic field. The plateau has a Hall conductivity quantized at ^f-^ with Co the 
c-axis lattice constant. We analyze the three-dimensional Hofstadter problem of a realistic tight- 
binding Hamiltonian for graphite, find the gaps in the spectrum, and estimate the critical value 
of the magnetic field above which the Hall plateau appears. When the Fermi level is in the bulk 
Landau gap, Hall transport occurs through the appearance of chiral surface states. We estimate 
the magnetic field necessary for the appearance of the three dimensional quantum Hall Effect to be 
15.4 T for electron carriers and 7.0 T for hole carriers. 

PACS numbers: 72.25.-b, 72.10.-d, 72.15. Gd 



Recent advances in the fabrication of single graphene 
sheets as well as the striking initial experiments on the 
relativistic quantum Hall effect in graphene 0, have 
generated intense interest in this remarkable material. 
Most of the theoretical and experimental research has fo- 
cused on the properties of the low energy excitations close 
to half filling which have a Dirac spectrum with a speed 
of light of the order of 10 6 m/s. The spin-unpolarized 
quantum Hall effect shows a sequence of plateaus at 
&xy = (n + \) x 4^- consistent with the existence of 
two Dirac cones as well as a spin degeneracy []], • 

The phenomenon of Hall conductivity quantization is 
however not restricted to two dimensions and can occur 
in bulk samples, albeit under more stringent conditions. 
It was first observed by Halperin [3] that for a three- 
dimensional (3D) electron system in a periodic potential, 
if the Fermi level lies inside an energy gap, the conduc- 
tivity tensor is necessarily of the form: 
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FIG. 1: (a) Graphite in Bernal stacking, (b) Under strong 
magnetic field, graphite is gapped in the bulk and exhibits 
chiral surface sheet states, (c) Idealized Brillouin zone for 
graphite, (d) Predicted 3D Hall conductivity, quantized in 
units of 1/cq. Only one plateau is observable in graphite. 
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where is the fully antisymmetric tensor and G is a 
reciprocal lattice vector of the potential which may or 
may not be zero. The periodic potential can be generated 
by either an underlying lattice or spontaneously by many- 
body effects through the formation of density waves. 

The existence of the three dimensional quantum Hall 
effect (3DQHE) has not, however, been observed in any 
bulk materials, although the effect has been observed in 
stacked superlattices of 2D semiconductor layers, each 
exhibiting a 2DQHE, where the inter-layer tunneling is 
engineered away [4| . The only current bulk material can- 
didates for 3DQHE are the Bechgaard salts [5|, The 
difficulty inherent in realizing the 3DQHE in a specific 
bulk material can be understood by considering the spec- 
trum of an electron in a periodic potential under an ap- 
plied magnetic field. In 2D, this Hofstadter problem [7[ 
results in magnetic energy bands. Associated with each 



band is an integer quantized Hall conductance; when the 
Fermi level lies in a gap, the observed Hall conductance 
is a sum over the contributions from each filled band. In 
3D, the momentum component parallel to the magnetic 
field disperses freely and the very existence of such gaps 
is not guaranteed; indeed it is exceptional. Tilting the 
magnetic field away from the crystallographic axes helps 
[8[, but does not guarantee, the opening of gaps. 

In this paper we show that a true, bulk 3DQHE is 
realized in graphite under a large magnetic field paral- 
lel to the c-axis. Three factors conspire to render this 
possible: the large Landau gap of the integer quantum 
Hall state in graphene, the weak interplane hopping, and 
the Bernal stacking. We first give a physical argument 
for the existence of 3DQHE based on adiabatic conti- 
nuity, then perform a full Hofstadter calculation in 3D 
of the band and surface states structure, using the realis- 
tic Johnson-Dresselhaus {<J [To[ Hamiltonian for graphite, 
plus a magnetic field. We find the minimum magnetic 
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FIG. 2: Spectrum of the zero mode and first few Landau 
levels of our toy model in B = 10 T interpolated between 
graphene (x = 0), where they do not k z disperse, and graphite 
(x = 1). The zero mode (blue) is doubly degenerate, giving 
a 3D Hall conductance of \— per independent K point per 
spin. 



field necessary for a 3DQHE to be 15.4 T for electrons 
and 7.0 T for holes, and show that only one Hall plateau 
(Fig. [T]) will be observed due to band closing in the higher 
Landau levels (LLs). Besides the obvious prediction of a 
plateau in off-diagonal conductive response, we also pre- 
dict that a correlated chiral surface state [ll|, occurs at 
the boundary of the sample. Based on the recent experi- 
mental focus on graphite [12[ we believe our prediction is 
testable with current experimental techniques. Previous 
Hall plateaus observed in graphite [HI, Q3] are either even 
in £?, or come in multiple sequences consistent with the 
graphene QHE plateaus [3, EEf , and are hence different 
than our prediction of a single quantized Hall plateau of 
Hall conductivity twice as large as the one observed in 
graphene for Fermi level in the first Landau gap. 

Graphite is a layered material consisting of weakly 
coupled parallel graphene layers in an ABABAB con- 
figuration known as Bernal stacking [16], as depicted 
in Fig. [TJ Let us initially set the inter-layer hopping 
to zero and place a large number of uncoupled stacked 
graphene layers in a strong magnetic field perpendicular 
to their surface. Each layer exhibits a relativist ic quan- 
tum Hall Effect as previously described with a LL energy 
E n = ±\ /r 2nhv/£, where £ = ^hc/eB is the magnetic 

length, n = 0, 1, and v « 10 6 m/s. Typical values of 

the gap are roughly 0.1 eV for B — 10 T and > 0.25 eV for 
B > 40 T, thus making graphene the first system to ex- 
hibit quantized Hall conductance at room-temperature. 
If we place the Fermi level in the first Landau gap, our 
system of uncoupled graphene layers trivially exhibits a 
3DQHE, with the bulk enveloped by a sheath of chiral 
surface states as in Fig. [TJ For uncoupled layers, there 



is no k z dispersion. We now adiabatically turn on intra- 
layer hopping, which causes almost all LLs to disperse 
with k z . The exception are the zero modes which are 
stable in an idealized model due to particle-hole sym- 
metry, as we show below. By adiabatic continuity, the 
sheath of chiral states (and hence the 3DQHE) must be 
stable as long as the Landau gap does not collapse. 

We first show that the statement above is true for 
a simplified toy- model of graphite. If we introduce 
a parameter x which interpolates between uncoupled 
graphene layers (x = 0) and graphite (x = 1), the 
minimal Hamiltonian of our system is a 4 x 4 Hermi- 
tian matrix acting on the lattice spinor (^ A , V^b > > ) > 
where (A,B) and (C,D) are the two inequivalent sites 
of the lower and upper graphene layers in the Bernal 
stacking, respectively. Its nonzero independent ele- 
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—2xt ± cos(t^Cq). The in-plane hopping is ty = 3.16 eV, 
t± = 0.39 eV is the inter layer (A-C) hopping in graphite 
(x = 1) [16], and a = 2.456 A is the graphene lat- 
tice constant (in-plane A- A distance). The c-axis lattice 
constant is c = 6.74 A (the distance between consec- 
utive lattice planes equal to |c ), and k± = k x ± ik y . 
The in-plane dispersion is expanded about the point 
K = (^,0,0), which is the Dirac point in graphene. 
In a uniform magnetic field, adopting the symmetric 
gauge A = ~B(—y,x : 0), the Kohn-Luttinger substitu- 
tion k± — > k x ± iky + j[^(A x ± iA y ) [17[. In the specific 
case of (x = 1) this model is the same as the one used in 
[Tit . Let b\ b be creation and annihilation Landau Level 
operators (with [6,6^] = 1), and introduce the notation 
c z = cos (^ k z c ). The Hamiltonian is 
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with e = B/B and B = (he/ e) /Snal = 7275 T, and is 
diagonalized in a basis i/j = \n + 1), |n), \n — 1)), for 
n > 0. We obtain below the 4 LL energies E n as 
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This spectrum has explicit particle-hole symmetry. For 

1/2 

n = there are eigenvalues at ± (e 2 t 2 +4x 2 t\ c 2 ) and 

a doubly degenerate level at E = 0. (All levels receive 
an additional double degeneracy owing to the existence 
of the inequivalent K point.) The next Landau bands 
are the lowest two energy levels of n = 1. Using this 
spectrum we find that the gap between the zero mode 
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and the LL above or below it cannot collapse upon in- 
terpolating between x = and x = 1 for any value of 
the magnetic field. As such, by adiabatic continuity, the 
Hall conductance when the Fermi level is in the 3D gap 
(with doubling for the two spins) is: 



TABLE I: Tight-binding parameters 
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For the Fermi level between the zero mode and the first 
LL, n = and a xy = ±4e 2 //ico- The Bernal stack- 
ing of graphite accounts for the extra factor of 2 rela- 
tive to graphene and J I m {^k^ I ^fc y ^) = 2n + 1 is 
the TKNN integer pjl of the relativistic graphene bands 
when the Fermi level is placed in the n th bulk gap. We ob- 
serve that the gap between the first and the second LL in 
graphite closes for any realistic value of the B field (Fig. 
[2]), and hence higher n plateaus will not be observed. Zee- 
man splitting could give rise to a a xy = plateau, but is 
smeared by the dispersion of the zero mode in the realis- 
tic graphite model used below, and hence the predicted 
Hall conductance is sketched in Fig. [TJ 

The existence of a full gap in the 3D LL spectrum is 
an artifact of the toy model involving only and t ± . It 
derives from the presence of a flat, twofold degenerate 
band at the Fermi level along the HK spines of the Bril- 
louin zone. The Bernal stacking is crucial, for if the indi- 
vidual graphene planes were stacked identically, all LLs 
would acquire a /^-dependent contribution to their en- 
ergies — 2t± cos(/c 2 c ), and would then overlap at all but 
extremely high magnetic fields [20]. More realistic mod- 
els for graphite, such as the Slonczewski-Weiss-McClure 
(SWMC) model [U or that of Johnson and Dressel- 
haus (JD) contain small c-axis B-B (D-D) hop- 

ping terms, through the open hexagons of the CD (AB) 
plane. Their value is small on the scale of nearest neigh- 
bor hopping - only 10 meV, leading to a bandwidth of 
40 meV along the HK spine. But the presence of such 
terms is crucial toward understanding the properties of 
graphite. At B = 0, they result in semimetallic behav- 
ior, whereas the toy model incorrectly predicts a zero 
gap semiconductor. For weak fields, they lead to over- 
lap of the LLs and destruction of the QHE. However, as 
we shall show, the principal gaps surrounding the central 
n = LLs survive for B > 7.0 T (holes) and B > 15.4 T 
(particles). Lightly doped graphite, then, will exhibit a 
3DQHE at these fields. We next describe our solution of 
the JD model in a magnetic field on a torus and a Laugh- 
lin cylinder, finding the LLs and the surface states, and 
determining the critical fields B c at which energy gaps 
open across the entire Brillouin zone [23]. 

The JD model 0, is a tight-binding Hamiltonian 
derived from the k • p theory of SWMC. Its nine param- 
eters are given in Tab. H In addition to nearest neigh- 
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FIG. 3: Definition of JD hopping matrix elements and flux 
assignment. Sites A and B belong to different graphene layers 
than sites C and D, and Bernal stacking corresponds to sites 
A and C differing by a c-axis translation. Not shown are the 
further- neighbor in-plane hoppings t^g and t^g. = 2np/q 
is the magnetic flux per hexagon in units of hc/e, and n runs 
from 1 to q, the number of hexagons in a magnetic unit cell. 



bor hoppings, there are also further neighbor hoppings, 
both in-plane (extending to third and fourth neighbor) 
and between planes. There is also an energy asymmetry 
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We then introduce a magnetic field 
through Peierls phases along the links, consistent with all 
rotational (and screw axis) symmetries of the underlying 
lattice, as depicted in Fig. 03 

We solve the model through a combination of exact di- 
agonalization, Lanczos method (for q > 1000, where the 
flux per hexagonal plaquette is 1/q Dirac quanta </> — 
hc/e), and low field expansion. For the bulk band struc- 
ture, we impose doubly periodic (i.e. toroidal) boundary 
conditions in the (x, y) plane, while to study edge (sur- 
face) states we impose singly periodic (i.e. cylindrical) 
boundary conditions. The latter case may also be solved 
via a transfer matrix formalism [24] . Due to the algebraic 
complexity of the problem, we postpone the explicit de- 
tails of the 5-field Hamiltonian and transfer matrix for a 
future publication. For q > 20 the Hofstadter broadening 
becomes negligible and the band energies as a function 
of k x become non-dispersive, corresponding to the real 
situation in which the magnetic field splitting is small 
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FIG. 4: Clockwise from upper left: (a) B = 5T, no gap 
present in the full spectrum, (b) B = 12 T, clear gap in the 
hole LL spectrum, (c) B = 20 T clear gaps for both hole and 
electron LL. Spin- up (blue) and spin-down (red) bands are 
shown, (d) Principal energy gaps surrounding n = LLs, 
including effects due to Zeeman splitting. The particle gap 
collapses at B e c = 15.4 T and the hole gap at = 7.0 T. All 
energies have been shifted upward by 100 meV. 



compared to the in-plane hopping amplitude. 

We observe the following characteristics of the spec- 
trum. The zero mode has a dispersion of 40 meV in the 
low field limit. The first LL, however, disperses strongly 
as a function of k z and for each value of q we scan the 
energy E(k z ) spectrum to look for the smallest gap. Fig. 
2] shows the bulk Landau Level spectrum for B = 5, 12, 
and 20 Tesla, as well as the field-dependence of the princi- 
pal gaps surrounding the central, weakly dispersive n = 
LLs. We find that both gaps are indirect for B < 30 T. In 
Fig. [5]we plot the bulk band and surface state spectrum 
as a function of k x and k z , which proves the existence 
of gaps and surface states over the entire Brillouin zone. 
On the Laughlin cylinder, for each value of k z , we ob- 
tain 2 edge states on each of the upper/lower edge of the 
cylinder. Unlike the low-field bulk LLs, the edge states 
disperse as a function of k x and cross the bulk gap to give 
2e 2 /h c per spin. We have checked that our numerical 
results match the theoretical analysis of the continuum 
k approximation of the JD model. 

At low fields, each LL accommodates ^a\h >/40 o = 
3.16 x 10 _6 £?[T] states per carbon atom. Accounting 
for the quadruple degeneracy of the LLs (two K points 
and two spin polarizations), the central n = levels will 
be filled for fields below 20 T at a doping of only 0.025%. 
For the lowest field for which we predict the effect, the 
doping is a modest 0.01%. Unlike in a many-body gap, 
the Fermi level is not pinned and the width of the LL 




FIG. 5: Zero mode, first electron Landau level and sur- 
face states spectrum over the full Brillouin zone for one spin 
species at B = 40 T. (a) Surface states spectrum (b) Zero 
mode and first electron Landau level spectrum (c) Surface 
states and bulk Landau levels together. The n — 1 (n = 0) 
levels are shifted by +0.1 eV (— 0.1 eV) for clarity. There are 
two surface states per spin species crossing the bulk gap. 



will be given by the width of the mobility single-particle 
gap in disordered graphite. Strong disorder leads to wide 
Hall plateaus, and weak disorder to narrow ones [2H ] . 

Besides the prediction of a quantized response, the chi- 
ral surface sheet should exhibit ballistic in-plane longitu- 
dinal response p xx — > as T — > 0. Remarkable transport 
properties are also present along the direction parallel to 
the magnetic field [111 ] . In this direction the system is a 
stable metal with a temperature independent resistivity 
p zz , with a value which can be much larger than h/e 2 , 
although, unlike previously predicted [11] we do expect 
the metallic phase to be unstable to very strong disor- 
der and impurity concentration which levitate the bulk 
extended states above the Fermi level (26| . 

In conclusion, we predict the appearance of one 
plateau of 3D quantized Hall transport in doped bulk 
graphite under strong magnetic fields parallel to its c- 
axis. We analyzed the Hofstadter problem in graphite 
and estimate the minimum field to be approximately 
7.0 T (for holes) and 15.4 T (for electrons). While our 
general conclusions are robust, differences among mate- 
rial parameters reported in the literature can shift the 
critical fields by several Tesla. The recent advances and 
interest in graphene and graphite make our prediction 
testable with current experimental techniques. 
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